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Recent interest in formation flying satellite systems has spurred a considerable amount of 
research in the relative navigation and control of satellites. Development in this area has 
included new estimation and control algorithms as well as sensor and actuator development 
specifically geared toward the relative control problem. This paper describes a simulation 
facility’, the Formation Flying Test Bed (FFTB) at NASA Goddard Space Flight Center, 
which allows engineers to test new algorithms for the formation flying problem with relevant 
GN&C hardware in a closed loop simulation. The FFTB currently supports the inclusion of 
GPS receiver hardware in the simulation loop. Support for satellite crosslink ranging 
technology is at a prototype stage. This closed-loop, hardware inclusive simulation capability 
permits testing of navigation and control software in the presence of the actual hardware 
with which the algorithms must interact. This capability provides the navigation or control 
developer with a perspective on how the algorithms perform as part of the closed-loop 
system. In this paper, the overall design and evolution of the FFTB are presented. Each 
component of the FFTB is then described. Interfaces between the components of the FFTB 
are shown and the interfaces to and between navigation and control software are described. 
Finally, an example of closed-loop formation control with GPS receivers in the loop is 
presented. 


1. Introduction 

S atellite formation flying systems have attracted much attention in recent years due to potential advantages over 
traditional, monolithic satellites. These advantages include increased effective aperture of on-orbit observing 
systems, capability to distinguish between spatial and temporal variations of in-situ measurements of the space 
environment, graceful degradation of the system as a result of spacecraft failures, system reconfiguration via a 
change in formation, phased formation population, and economies of scale in the construction of spacecraft. In fact, 
the President’s Commission on Implementation of United States Space Exploration Policy recently identified 
formation flying as one of seventeen enabling technologies for attainment of exploration objectives [1], 

NASA has several formation flying systems in its plan for upcoming missions including the Micro-Arcsecond X- 
ray Imaging Mission [2], Stellar Imager [3], Terrestrial Planet Finder [4], Submillimeter Probe of the Evolution of 
Cosmic Structure [5], Constellation-X [6], and Magnetospheric MultiScale Mission [7], Additionally, one of the 
five concepts that the New Millennium Program is considering for the technology’ demonstration mission ST-9 [8] is 
a precision formation flying mission. 

Among the most challenging aspects of these missions is the relative navigation and control of spacecraft within 
a formation. Many researchers have developed novel concepts and algorithms for formation navigation and control 
over the past few years. The aim of the FFTB, shown in Figure 1, is to provide an environment where these new 
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algorithms can be tested in the presence of real hardware and where hardware can be tested against the unique 
challenges of formation flying requirements. Major benefits of this type of testing include the exercise of hardware- 
software interfaces, early identification of hardware problems, and realistic assessment of algorithm performance in 
the presence of actual hardware error characteristics. 



Figure 1. The Formation Flying Test Bed 


II. FFTB Architecture 

The FFTB hardware architecture [9] is shown in Figure 2. The simulation originates at the Environment 
Computer, which generates vehicle trajectories via high fidelity orbit propagation for each spacecraft in formation. 
The trajectories generated by the Environment are used as truth for the simulation. The Environment Computer then 
supplies ep'nemeris information to the Spirent system in real-time by monitoring a one pulse-per-second (PPS) signal 
originating from the GPS Simulator as shown or an external source such as a cesium oscillator The GPS Simulator 
uses that information together with almanac data from the GPS constellation and the receivers’ antenna 
configurations to produce up to four Radio Frequency (RF) signals representative of the aggregate signals from the 
visible portion of the GPS constellation (up to 16 GPS satellites per RF signal) for each GPS receiver in the 
formation. Those RF signals are routed via coaxial cable into the GPS receiver hardware where measurements are 
generated and sent via RS-232 interface (serial) to the Flight Computers. The Flight Computers host the navigation 
and control processes under the management of a top level process called the Flight Executive that manages 
incoming data, calls the navigation process, and uses the navigation results as input to the formation control process. 
The Flight Executive then feeds the maneuver commands from the control process back to the Environment 
Computer for incorporation into the equations of motion, thus closing the control loop. The Crosslink Channel 
Simulator (CCS), a current development at the prototype stage [10], is also shown in Figure 2. As shown, the CCS 
will enable insertion of satellite crosslinks in the simulation loop. The CCS is the analog to the GPS Simulator for 
crosslink hardware. That is, on the basis of spacecraft trajectory information received from the Environment 
Computer, it modifies inter-spacecraft RF transmissions, introducing appropriate levels of time delay, Doppler 
frequency shift, and signal attenuation, thus enabling realistic simulation of inter-spacecraft communication and 
ranging. Both the Flight Computer and the Environment supply information to the Visualization Computer for data 
plotting and near real-time, three-dimensional animation of the formation through Satellite Tool Kit (STK™). 

III. FFTB Elements 


A. Environment Computer 

The Environment Computer is a standard personal computer running the Linux operating system. A kernel 
device driver listens to the incoming PPS in order to synchronize the simulation with the GPS Simulator. The 
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spacecraft orbits are propagated in real-time using the Spacecraft Trajectory and Attitude Real-time Simulation 
(STARS). STARS currently uses the orbit propagation capability resident in the GPS Enhanced Onboard Navigation 
System (GEONS) [11], and the attitude profile is prescribed as Earth-pomtmg. State updates are sent to the GPS 



Simulator at 10 Hz; consequently, the orbit propagation step size is 0.1 s. In the future the step size will be 
increased in order to avoid long-term numerical error accumulation, and intermediate state updates will be 
calculated using an interpolation scheme. STARS has the capability to accept control maneuvers over UDP 
(Universal Datagram Protocol) Socket connection from the Flight Computer. The maneuver messages are currently 
sent as AV (instantaneous change in velocity) values in the Radial, 

Transverse, Orbit-Normal (RTN) frame of the spacecraft where radial is 
the unit vector of the position vector, transverse is perpendicular to the 
position vector in the orbit plane and in the same sense as the velocity 
vector, and orbit normal completes the right handed triad. The AV 
message is converted to a constant acceleration applied over the 
integration step. 

B. GPS Simulator 

The GPS Simulator is a Spirent STR4760 (Figure 3) with four 
available RF outputs which can accommodate up to four GPS receivers 
with one antenna each. The antenna pattern for each receiver is 
configurable; normal FFTB procedure is to assign hemispherical 
coverage. The GPS Simulator can simulate several sources of GPS 
measurement error including ionospheric delays, multipath, and Selective 
Availability. Up to sixteen GPS satellite signals are generated 
simultaneously on each RF output; default selection of the simulated GPS 
signals is by elevation angle but can be configured to use signal strength Figure 3. Spirent 4760 GPS 
which is useful for spacecraft in orbits which go above the GPS Simulator 

constellation. 

C. GPS Receivers 

The FFTB currently has an inventory of several space-capable GPS receivers. Four GSFC-developed PiVoT 
receivers [12] are available and are currently being updated for above-the-constellation use. Two Orion receivers 
from the German Space Operations Center are also available and include a useful relative navigation function 
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whereby they share GPS measurements via direct RS-232 connection and the relative satellite state is estimated. A 
commercial off-the-shelf receiver, the Ashtech G-12 [13], is also available for performance baselining. 


D. Flight Computers 

The Flight Computers host the navigation 
software process and formation controller as showm 
in Figure 4. The formation controller is a modular 
component that can accommodate either 
MATLAB™ scripts or standard C/C++ code. As 
mentioned previously, the crosslink interface is not 
yet implemented, but, with the delivery of the 
operational CCS, it is expected to be in place before 
February 2005. In lieu of crosslink hardware in the 
simulation, the Flight Executive is capable of 
passing inter-spacecraft data messages to the other 
Flight Computers via UDP sockets. The Flight 
Executive is a process that sits on top of the 
navigation and control functions and is responsible 
for data management and execution of the 
navigation and control functions. 

E. Visualization Computer 

The Visualization Computer receives the truth data 
from the Environment Computer and navigation 
results from the Flight Computers, both via UDP 
socket connection. Truth data are then passed to 
STK™ for animation (Figure 5) of the formation. The 
Visualization Computer also has a near-real-time 
plotting capability which provides easy access to 
important information like the number of GPS 
satellites tracked, navigation errors, etc. 
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Figure 4. Flight Executive Schematic 



F. Crosslink Channel Simulator 

The CCS will enable FFTB users to insert Figure 5. STK Visualization 

crosslink hardware in the simulation for more realistic 
inter-spacecraft communication and ranging. It 
currently exists in the prototype state as shown in 
Figure 6. The capability to introduce signal delay, 

Doppler shift and attenuation onto a crosslink signal 
has been demonstrated. Delivery of the operational 
version of the CCS is expected before February 2005. 

The operational version of the CCS will introduce 
signal delay with resolution down to 30 ps at a 
maximum range of 64,000 km. 

IV. Formation Flying Test Bed Simulator 

An all-software, non-realtime simulation has been 
developed to mimic the performance of FFTB real-time, 
hardware-in-the-loop simulation. This simulation, 
henceforth refered to as FFTB Sim, is a software version 
of the FFTB which provides an interface identical to the Flightexec-Matlab™ interface to allow users to test FFTB 
guidance and control code in non-real-time. The simulation is driven by the FreeFlyer™ orbit software, and 
includes high-fidelity gravity, drag, and solar radiation pressure dynamics, as well as measurement noise affecting 
the estimated states. 



Figure 6. CCS Prototype 
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V. Example Result 

To demonstrate closed loop control for a formation flying mission, we perform a simulation in the FFTB that 
includes an actively controlled spacecraft and a passive spacecraft, each equipped with an Orion GPS receiver. 
Absolute orbit determination is performed on the active spacecraft flight computer using the GPS Enhanced 
Onboard Navigation System (GEONS), which processes pseudorange data from one of the Orion GPS receivers in 
an Extended Kalman Filter (EKF). For this scenario, GEONS estimates the receiver clock error bias and drift, 
absolute position and velocity, and drag coefficient. Control accelerations are handled in the EKF by including the 
accelerations in the state propagation, and by increasing the position and velocity covariance whenever control is 
applied. 

Relative navigation is performed by the Orion receivers, which exchange raw measurements over a serial port, and 
output time and relative position and velocity with respect to the RTN frame of the host satellite. The relative 
navigation algorithm is described in more detail in [14] and results with a similar configuration of the Orion receiver 
can be seen in [15]. 

Control is computed for the active spacecraft using a Matlab™ function call commanded by the Flight Executive, 
which has, as inputs, time, absolute and relative states, and, as outputs, spacecraft control acceleration in the Radial, 
In-Track, Cross-Track (RIC) frame**. 

The control law implemented in this study is a simple proportional-derivative feedback of the difference between the 
desired trajectory and the relative state information provided by the relative navigation system. Absolute state 
estimates are not used directly in the control feedback, but are available for computation of coordinate 
transformations. Thrust is assumed to be equally available in any direction, with magnitude limited to lOmN. 
Control is calculated and applied at a frequency of 1Hz. The desired trajectory is circular motion of one spacecraft 
about the other with a fixed separation, and a time varying phase angle within the plane of relative motion of the 
circular formation. 

The two spacecraft formation is injected into a nearly circular parking orbit with mean semi-major axis of 6823 km, 
and mean inclination of 28 degrees. The two spacecraft are separated by a 1km in-plane separation in the parking 
orbit, with the maneuvering spacecraft trailing. The transfer to the initial precision formation flying configuration is 
performed using two maneuvers separated by half the orbital period. The first maneuver is a combination radial and 
out-of-plane bum, which puts the maneuvering spacecraft on course to a point 100 meters ahead of and 100 meters 
out-of-plane from the passive spacecraft. The second maneuver is a radial bum half an orbital period after the first 
bum to achieve a natural circular formation. Closed-loop control of the precise circular motion trajectory is 
simulated in FFTBSim with relative position and velocity noise of lm and 0.5 cm/s, and in the FFTB hardware-in- 
the-loop (HWIL) simulation described above. 


Table 1. Hardware-in-the-loop steady-state estimation errors for a 4 hour simulation 


Navigation Type 

Position Estimation Error [m] 

Velocity Estimation Error [cm/s] 


Mean 

Max 

Mean 

Max 

GEONS absolute state estimate 

1.08 

2.36 

0.218 

0.494 

Orion abs. state est. (maneuvering s/c) 

1.29 

8.67 

3.058 

22.162 

Orion abs. State est (passive s/c) 

0.23 

35.93 

2.772 

40.503 

Orion relative state estimate 

0.99 

8.06 

0.510 

2.400 


Absolute and relative state estimation errors from a 4-hour hardware-in-the-loop FFTB simulation of closed loop 
formation control are presented in Table 1. The table presents steady-state values from the final 3 hours of the 
simulation. The GEONS absolute state estimation results are as expected, with position errors on the order of a 
meter, and velocity errors on the order of tenths of cm/s. The relative state error of the Orion receivers is 


All relative velocities described in this paper are assumed to be with respect to the rotating reference frame, and 
therefore include the “omega-cross” term. 
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considerably larger than expected. Comparing absolute state errors from the two Orion solutions, we see that the 
state error is considerably larger for the maneuvering spacecraft than for the passive spacecraft. In hindsight this 
result is not surprising. We are seeing the effect of noisy control application on the navigation accuracy. The 
receiver has no knowledge of the control being applied, and cannot be expected to perform as well in such a 
perturbed environment. There would be a similar effect in the GEONS output accuracy if the filter had no 
knowledge of the control acceleration being applied. 

Table 2 shows controlled inter-spacecraft range errors from both all-software (FFTBSim) and FFTB hardware-in- 
the-loop (HWIL) simulations. Initial condition errors in this simulation are due to navigation and thrust performance 
errors in the transfer from the parking orbit to the desired formation configuration. 

Figure 7 compares the control performance for the software-only and hardware-in-the-loop simulations. The 
HWIL mean range error is about 2.5 times worse than the FFTBSim result. This inconsistency is due to a number of 
effects, including modeling of the state estimation error as white Gaussian noise in FFTBSim, and inconsistent real- 
time performance in the Flight Executive. 

Perhaps the most notable data from Table 2 are the I20j . . — ■ — 

thruster duty cycles and total AV values. Thrust is being . _ controlled (Software) 

applied about half as often in the HWIL simulation as in 115 — Controlled (HWIL) 

the software-only simulation. Since the Flight 

Executive software is not synchronized, the control law 110 • \ 

is not being executed at exactly once per second, ]=~ / \ / 

resulting in impaired performance. ^ j \ ! \ 

VI. Conclusions / \ / * . 

The FFTB is a viable and useful tool for assessing \ _/ *, 

the efficacy of formation navigation and control _ \ 

algorithms in the presence of real navigation hardware. 9b ' ' 

Moreover, by effectively exercising hardware-software 

interfaces and exposing algorithms to real hardware 9 Cp ^ 

error characteristics, the FFTB serves as the "relevant (jme j s j 

environment" for Technology Readiness Level 

advancement. Fioure 7. Controlled and uncontrolled inter- 


time [s] 

Figure 7. Controlled and uncontrolled inter- 
spacecraft range, with proportional derivative 
control targeting 100m separation 


Table 2 Control performance results 


Simulation Type 



Inter-Spacecraft Range Error [cm] 

Total AV 

Mean 

Max 

[m/s] 

13.24 

44.92 

1.396 

31.14 

126.30 

0.777 


Thruster Duty Cycle 
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